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We review how RHIC is expected to deepen our understanding of the spin structure 
of longitudinally and transversely polarized nucleons. After briefly outlining the 
current status of spin-dependent parton densities and pointing out open questions, 
we focus on theoretical calculations and predictions relevant for the RHIC spin 
program. Estimates of the expected statistical accuracy for such measurements 
are presented, taking into account the acceptance of the RHIC detectors. 


1. Lessons from (Un)polarized DIS 

Before reviewing the prospects for spin physics at the BNL-RHIC we briefly 
turn to longitudinally polarized deep-inelastic scattering (DIS) and what we 
have learned from twenty years of beautiful data 1 . Figure 1 compares the 
available information on the DIS structure function gi(x,Q 2 ) to results of a 
typical next-to-leading order (NLO) QCD fit. From such types of analyses a 
pretty good knowledge of certain combinations of different quark flavors has 
emerged, and it became clear that quarks contribute only a small fraction to the 
proton’s spin. However, there is still considerable lack of knowledge regarding 
the polarized gluon density A < 7 , which is basically unconstrained by present 
data, the separation of quark and antiquark densities and of different flavors, 
and the orbital angular momentum of quarks and gluons inside a nucleon. In 
addition, spin effects with transverse polarization at the leading-twist level, 
the so-called ‘transversity’ densities, have not been measured at all. With 
the exception of orbital angular momentum RHIC can address all of these 
questions as will be demonstrated in the following 2 . 

There is also an important difficulty when analyzing polarized DIS data in 
terms of spin-dependent parton densities: compared to the unpolarized case 
the presently available kinematical coverage in x and Q 2 and the statistical 
precision of polarized DIS data are much more limited 1 . As a consequence, 
one is forced to include data into the fits from (x 7 Q 2 )-regions where fits of 
unpolarized leading-twist parton densities start to break down, see Fig. 1. 
Data from RHIC, taken at ‘resolution’ scales Q 2 where perturbative QCD and 
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Figure 1. Available information on g\ (x. Q' 2 ) as collected by fixed-target experiments 1 
compared to results of a NLO QCD fit (solid lines). The indicated rectangular and tri¬ 
angular regions contain data which would not pass kinematical cuts of Q 2 > 4 GeV 2 and 
W 2 > 10GeV 2 , respectively, typically imposed in all fits to unpolarized DIS data. 

the leading-twist approximation are supposed to work, can shed light on the 
possible size of unwanted higher-twist contributions in presently available sets 
of polarized parton distributions. 

2. Spin Physics at RHIC with Longitudinal Polarization 
2.1. Prerequisites 

The QCD-improved parton model has been successfully applied to many high 
energy scattering processes. The predictive power of perturbative QCD follows 
from the universality of the parton distributions and fragmentation functions. 
Once extracted from data they can be used to make definite predictions for 
other processes. This property is based on the factorization theorem where a 
cross section is written as a convolution of perturbatively calculable partonic 
hard scattering coefficients da^ b and appropriate parton densities f a ^ and/or 
fragmentation functions . To be specific, let us consider the inclusive pro- 
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duction of a hadron H , e.g., a pion, in unpolarized proton-proton collisions: 

dr 


d(j^ f dd^ 

= dx a dx b dzf a (x a ,p)f b (x b ,p)—^-(x a ,x b ,z,r,p)D^(z,p). (1) 


dT 


abc 


Here, T stands for any appropriate set of kinematical variables like the trans¬ 
verse momentum pr and/or rapidity y of the observed hadron. The functions 
f a , b and D^ embody non-perturbative physics. However, once they are known 
at some initial scale po, their scale /z-dependence is calculable perturbatively 
via a set of evolution equations. The factorization scale /x, introduced on the 
r.h.s. of (1), separates long- and short-distance phenomena, p is completely 
arbitrary but usually chosen to be of the order of the scale characterizing the 
hard interaction, for instance pt■ Since the l.h.s. of (1) has to be independent 
of p (and other theoretical conventions), any residual dependence of the r.h.s. 
on the actual choice of p gives an indication of how well the theoretical calcu¬ 
lation is under control and can be trusted. In particular, leading order (LO) 
estimates suffer from a strong, uncontrollable scale dependence and hence are 
not sufficient for comparing theory with data. Figure 2 shows a typical factor¬ 
ization scale dependence for various processes and experiments as a function of 
Pt- Clearly, the situation is only acceptable at collider experiments where one 
can access pt ^ 5 GeV <C y/S /2 = p^ ax (\/S is the available c.m.s. energy), pt 
values of about 1-2 GeV, accessible at fixed-target experiments, are not suffi¬ 
cient to provide a large enough hard scale to perform perturbative calculations 
reliably. For simplicity we have not distinguished between renormalization and 
initial/final-state factorization scales in (1) which can be chosen differently. 



Figure 2. Typical factorization scale dependence for various processes and experiments as 
a function of pT- Shown is the cross section ratio for two choices of scale, pt and Pt/2- 
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So far we have neglected the spin information contained in parton densities 
(and fragmentation functions). Eq. (1) can be easily extended to incorporate 
polarization by replacing all unpolarized quantities by their spin-dependent 
counterparts, like, for instance, f a ^ —> A f a ^ and dd c ab /dT —» dAd^ b /dT. If a 
hard-scattering process with incoming protons having definite spin orientations 
is studied, as at RHIC, one gains access to the spin distributions of quarks and 
gluons in a longitudinally (or transversely) polarized proton. In practice, spin 
experiments measure not the polarized cross section, dAa/dT, itself, but the 
spin asymmetry, which is given by the ratio of the polarized and unpolarized 
cross sections, e.g., for our example, Eq. (1), it reads 


H _ dAa H /dT 
LL = d<r H /dT 


( 2 ) 


To denote the type of polarization of the colliding hadrons in (2) we use the 
subscripts ‘L=longitudinal’ and ‘T=transverse’. At RHIC one can also study 
doubly transverse spin asymmetries, Att, and single spin asymmetries Al, At 
( the latter is often called An) where only one of the protons is polarized. 


2.2. Accessing A g 

The main thrust of the RHIC spin program 2 is to pin down the so far elusive 
gluon helicity distributions A g(x,g). The strength of RHIC is the possibility 
to probe Ag(x. ji) in a variety of hard processes 2 , in each case at sufficiently 
large pt where perturbative QCD is expected to work. This not only allows 
to determine the ai-shape of Ag(x. //) for x > 0.01 but also verifies the univer¬ 
sality property of polarized parton densities for the first time. In the following 
we review the status of theoretical calculations for processes sensitive to A g, 
experimental aspects can be found, e.g., in 3 . 

The ‘classical’ tool for determining the gluon density is high-px prompt 
photon production due to the dominance of the LO Compton process, qg —> 'yq. 
Exploiting this feature, both RHIC experiments, PHENIX and STAR, intend 
to use this process for a measurement of Ag. Apart from ‘direct’ mechanisms 
like qg —> 'yq, the photon can also be produced by a parton, scattered or created 
in a hard QCD reaction, which fragments into the photon. Such a contribution 
naturally arises in a QCD calculation from the necessity of factorizing final- 
state collinear singularities into a photon fragmentation function Df. However, 
since photons produced through fragmentation are always accompanied by 
hadronic debris, an ‘isolation cut’ imposed on the photon signal in experiment, 
e.g., a ‘cone’, strongly reduces such contributions to the cross section. 

The NLO QCD corrections to the direct (non-fragmentation) processes have 
been calculated in 4 and lead to a much reduced factorization scale dependence 
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Figure 3. All for prompt photon production in NLO QCD as a function of />■/■ for different 
sets of parton densities. The ‘error bars’ indicate the expected statistical accuracy for the 
PHENIX experiment. Figure taken from®. 

as compared to LO estimates. In addition, Monte Carlo codes have been 
developed 5,6 , which allow to include various isolation criteria and to study 
photon-plus-jet observables. The latter are relevant for A g measurements 
at STAR 2,3 . Since present comparisons between experiment and theory are 
not fully satisfactory in the unpolarized case, in particular in the fixed-target 
regime, considerable efforts have been made to push calculations beyond the 
NLO of QCD by including resummations of large logarithms 7 . It is hence 
not unlikely that a better understanding of prompt photon production can be 
achieved soon. Figure 3 shows A^ L as predicted by a NLO QCD calculation 6 as 
a function of the photon’s transverse momentum. The rapidity cut \r/\ < 0.35 
matches the acceptance of the PHENIX detector. The important result is that 
the expected statistical errors are considerably smaller than the changes in 
d ue 1° different spin-dependent gluon densities over a wide range of px- 
RHIC should be able to probe A g in prompt photon production. 

Jets are another key-process to pin down A g at RHIC: they are copiously 
produced at \/S = 500 GeV, even at high pr, 15 < pr < 50 GeV, and gluon- 
induced gg and qg processes are expected to dominate in accessible kinematical 
regimes. Due to limitations in the angular coverage, jet studies will be per¬ 
formed by STAR only. PHENIX can alternatively look for high-p-T leading 
hadrons, such as pions, whose production proceeds through the same partonic 
subprocesses as jet production. Hadrons have the advantage that they can be 
studied also at \f,S = 200 GeV and down to lower values in pt than jets as 
they do not require the observation of clearly structured ‘clusters’ of particles 
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Figure 4. As in Fig. 3 but now for tiigli-pv jet production. The ‘error bars’ are for the 
STAR experiment taking into account its acceptance. Figure taken from 8 . 

(jets). On the downside, they require a fragmentation function in the theoret¬ 
ical description, cf. Eq. (1), which is, however, fairly well constrained by e + e - 
data. It should be emphasized that in the unpolarized case, the comparison 
between NLO theory predictions with jet production data from the Tevatron 
is extremely successful. 

The NLO QCD corrections to polarized jet production are available as a 
Monte Carlo code 8 . Apart from a significant reduction of the scale dependence, 
they are also mandatory for realistically matching the procedures used in ex¬ 
periment in order to group final-state particles into jets. For single-inclusive 
high-pr hadron production the task of computing the NLO corrections has 
been completed only very recently 9 . Figure 4 shows All for single-inclusive 
jet production at the NLO level as a function of the jet pt- A cut in rapidity, 
|??| < 1, has been applied in order to match the acceptance of STAR. The 
asymmetries turn out to be smaller than for prompt photon production, but 
thanks to the much higher statistics one can again easily distinguish between 
different spin-dependent gluon densities. Very similar results are obtained for 
single-inclusive pion production 9 . 

The last process which exhibits a strong sensitivity to A g is heavy flavor 
production. Here, the LO gluon-gluon fusion mechanism, gg —> QQ, dominates 
unless pt becomes rather large. Unpolarized calculations have revealed that 
NLO QCD corrections are mandatory for a meaningful quantitative analysis. 
In the polarized case they have been completed only very recently for single- 
inclusive heavy quark production 10 . Again, one observes a strongly reduced 
scale dependence for charm and bottom production at RHIC energies. It turns 
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Figure 5. NLO single-inclusive charm production asymmetry (rescaled by l/a:™ m ) as a 
function of a;™ m = p™ m /p™ ax for different sets of parton densities. The ‘error bars’ are 
for the PHENIX experiment and include a detection efficiency for the channel c —> eX as 
modeled by PYTHIA. Figure taken from 10 . 

out that the major theoretical uncertainty stems from the unknown precise 
values for the heavy quark masses 10 . Since the heavy quark mass already sets 
a large scale, one can perform calculations for small transverse momenta or 
even for total cross sections which give access to the gluon density at smaller 
;c-values than relevant for jet or prompt photon production. 

Heavy flavors are not observed directly at RHIC but only through their 
decay products. Possible signatures for charm/bottom quarks at PHENIX are 
inclusive-muon or electron tags or ^e-coincidences. The latter provide a much 
better c/ 6 -separation which is an experimental problem. In addition, lepton 
detection at PHENIX is limited to \y\ < 0.35 and 1.2 < \y\ < 2.4 for electrons 
and muons, respectively. Since heavy quark decays to leptons proceed through 
different channels and have multi-body kinematics, it is a non-trivial task to 
relate, e.g., experimentally observed pr-distributions of decay muons to the 
calculated pr-spectrum of the produced heavy quark. One possibility is to 
model the decay with the help of standard event generators like PYTHIA 11 by 
computing probabilities that a heavy quark with a certain (pt,p) is actually 
seen within the PHENIX acceptance for a given decay mode. Figure 5 shows 
a prediction for the charm production asymmetry All at PHENIX in NLO 
QCD for the inclusive-electron tag. The sensitivity to A g is less pronounced 
than for the processes discussed above. It remains to be checked if heavy 
flavor production at RHIC can be used to extend the measurement of A g 
towards smaller x-values. Also, progress has to be made to solve the long¬ 
standing puzzle that the inclusive 6 -rate as predicted by QCD is too small in 
unpolarized pp, ep, and 77 collisions 12 . 















2.3. Further Information on A q and A q 

Inclusive DIS data only provide information on the sum of quarks and anti¬ 
quarks for each flavor, i.e., A q + A q. At RHIC a separation of Ait, Art, Ad, 
and Ad can be achieved by studying IR ± -boson production. Exploiting the 
parity-violating properties of VE ± -bosons, it is sufficient to measure a single 
spin asymmetry, A ^ , with only one of the colliding protons being longitudi¬ 
nally polarized. The idea is to study A^ as a function of the rapidity of the W, 
yw, relative to the polarized proton 13 . In LO it is then easy to show 13,2 that 
for VP + -production, ud —> W + , and large and positive (negative) yw, A™ is 
sensitive to A u/u (A d/d). Similarly, IT - -production probes A d/d and A u/u. 
The NLO QCD corrections for Al as well as the factorization scale dependence 
are small 14 . Experimental complications 2 arise, however, from the fact that 
neither PHENIX nor STAR are hermetic, which considerably complicates the 
reconstruction of yw- The anticipated sensitivity of PHENIX on the flavor 
decomposed quark and antiquark densities is illustrated in Fig. 6. 
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Figure 6. Expected statistical accuracy for A q/q from A l overlayed on two sets of parton 
densities. The full [open] circles refer to Al(W+) [Al(W~ )]. Figure taken from 2 . 

Semi-inclusive DIS measurements, ep —> HX, are another probe to sep¬ 
arate quark and antiquark densities. HERMES has recently published first 
preliminary results 15 . The accessible x-range for the A q and A q densities is 
comparable to that of RHIC, see Fig. 6, but at scales Q ~ 1 — 2GeV rather 
than Mw ■ The combination of both measurements can provide an important 
test of the QCD scale evolution for polarized parton densities. 
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2.4. Towards a Global Analysis of Upcoming Data 

Having available at some point in the near future data on various different 
reactions, one needs to tackle the question of how to set up a ‘global QCD 
analysis’ for spin-dependent parton densities. The strategy is in principle clear 
from the unpolarized case: an ansatz for the densities at some initial scale /z 0 , 
given in terms of some functional form with a set of free parameters, is evolved 
to a scale p relevant for a certain data point. A y 2 -value is assigned that 
represents the quality of the comparison of the theoretical calculation to the 
experimental point. The parameters are varied until eventually a global mini¬ 
mum in is reached mutually for all data points. In practice, this approach 
is not fully viable since the numerical evaluations of the cross sections in NLO 
QCD are usually time-consuming as they require several tedious integrations. 
Hence the computing time for a QCD fit easily becomes excessive. 

In the unpolarized case, the wealth of DIS data already provides a pretty 
good knowledge of the parton densities, and reasonable approximations can 
be made for the most time-consuming processes. For instance, one can ab¬ 
sorb all NLO corrections into some pre-calculated ‘correction factors’ AT, and 
simply multiply them in each step of the fit to the LO approximation for the 
cross sections which can be evaluated much faster. In the polarized case, it 
is in general not at all clear whether such a strategy will work. Here, parton 
densities are known with much less accuracy so far. It is therefore not pos¬ 
sible to use ACfactors reliably. In addition, spin-dependent parton densities 
as well as partonic cross sections may oscillate, i.e., have zeros, in the kine- 
matical regions of interest such that predictions at LO and the NLO can show 
marked differences. Clearly, in the polarized case the goal must be to find a 
way of implementing efficiently, and without approximations, the exact NLO 
expressions for all relevant hadronic cross sections. A very simple and straight¬ 
forward method based on ‘double Mellin transformations’ was proposed in 16 . 
Recently, its actual practicability and usefulness in a global QCD analysis has 
been demonstrated 1 ' in a case study based on fictitious prompt photon data. 

3. Spin Physics at RHIC with Transverse Polarization 

At RHIC one can also study collisions of transversely polarized protons 2 giving 
access to the completely unmeasured leading-twist ‘transversity’ densities 18 5f. 
Upon expressing transversely polarized eigenstates as superpositions of helicity 
eigenstates, 5f reveals its helicity-flip, chiral-odd nature which explains its 
elusiveness. Other striking features of Sf are that no transversity gluon density 
is possible for spin-1/2 targets and the fact that 6f(x,p) ‘evolves away’ at all 
x with increasing scale p. 
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3.1. Double Spin Asymmetries Att 

The requirement of helicity conservation in hard scattering processes implies 
that chirality has to be flipped twice in order to be sensitive to transversity. 
One possibility is to have both colliding protons transversely polarized and 
to study double spin asymmetries Att- Since gluons play an important or 
even dominant role in almost all unpolarized production processes, Att is 
expected to be very small in general 18 . In addition, Att is further dimin¬ 
ished by the requirement of a double chirality flip, which excludes some of 
the ‘standard’ 2 —> 2 amplitudes to contribute, whereas the remaining ones 
are color-suppressed. In principle, the most favorable reaction for determin¬ 
ing transversity is the Drell-Yan process, pp —> y + p~, which has no gluonic 
contribution in LO. However, a recent NLO study of upper bounds for Att, 
due to Soffer’s inequality 19 , 2|<5/| < / + A/, has revealed that the limited 
muon acceptance for PHENIX threatens to make a measurement of transver¬ 
sity elusive in this channel 20 . Recently is has been shown that, although Att 
is rather minuscule, jet and prompt photon production can be a useful tool to 
decipher transversity at RHIC 21 , see Fig. 7. It should be noted that NLO QCD 
corrections for these processes are still lacking. Needless to say that such mea¬ 
surements are challenging, albeit not completely impossible. The experimental 
finding of a much larger Att than theoretically expected would constitute a 
new spin puzzle. 
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Figure 7. Left: upper bound for Att for single-inclusive jet production at RHIC as a 
function of pt- Jot rapidities are integrated over the detector acceptance (—1 < 77 < 2 ). 
The shaded bands represent the uncertainties due to variations of the factorization scale in 
the range Pt/2 < p < 2 pj*. The expected statistical accuracy is indicated as ‘error bars’. 
Right: same as on the l.h.s. but now for prompt photon production. The photon rapidity 
has been integrated over the range |? 7 | < 0.35. Figure taken from 21 . 
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Alternative observables have only one transversely polarized initial hadron 
and a fragmentation process in the final state that is sensitive to transverse 
polarization 18 . Several processes have been identified as being potentially suit¬ 
able for a measurement of Sf: the production of transversely polarized A 
hyperons , the asymmetry in the py-distribution of a hadron in a jet around 
the jet axis 23 (‘Collins effect’), or the interference between s- and p -waves of a 
two pion system 24 (‘interference fragmentation’). Such measurements can be 
also carried out at HERMES and COMPASS. However, in all cases involving 
these novel fragmentation effects, the analyzing power is a priori unknown and 
may well be small, and often there are competing mechanisms for generating 
the same observable that do not involve transversity. 


3.2. Single Spin Asymmetries A 
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Figure 8. (a): for pion production in pp collisions at \/S = 20GeV compared to E704 

data. Predictions for RHIC 25 , VS = 200 GeV, for pt = 4 GeV are superimposed, (b): 
py-dependence of .4 n for RHIC at xp = 0.4. Figure taken from 2 . 

Surprisingly large single transverse spin asymmetries An have been ob¬ 
served over many years in low-energy fixed-target experiments 18 , e.g, for 
pp —> 7 tX. At the leading-twist level for standard (fep-integrated) parton 
distributions Ajy is exactly zero. One possible explanation is based on a gen¬ 
eralized factorization theorem in perturbative QCD, where a non-vanishing 
An can arise, for instance, as a convolution of some calculable hard-scattering 
function with an ordinary twist-two parton density from the unpolarized pro¬ 
ton and a new twist-three quark-gluon correlation function characterizing the 
polarized hadron 25 . In a simple model for these correlation functions, which 
are believed to give the dominant contribution to An in this approach, a qual¬ 
itative description of the available data is possible 25 , see Fig. 8, and various 
predictions have been made which can be tested at RHIC. In particular, at 
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RHIC one should see the fall-off of An with pt associated with its higher-twist 
nature, see Fig. 8. 

An alternative approach to An introduces intrinsic transverse momentum 
hr into distribution and fragmentation functions 18 . This opens a Pandora’s 
box with many new and unknown functions. Each of three possible mechanisms 
on its own can account for An data 18 . Needless to say that it is very difficult 
to disentangle all these effects. RHIC, with the help of other experiments, can 
help to shed some light on the origin of transverse single spin asymmetries. 

4. Exploring Physics Beyond the Standard Model 

Spin observables are also an interesting tool to uncover important new physics. 
One idea is to study single spin asymmetries A^ for large-pT jets. In the stan¬ 
dard model Al can be only non-zero for parity-violating interactions, i.e., 
QCD-electroweak interference contributions, which are fairly small. The exis¬ 
tence of new parity-violating interactions could lead to sizable modifications 26 
of A l . Possible candidates are new quark-quark contact interactions, charac¬ 
terized by a compositeness scale A. RHIC is surprisingly sensitive to quark 
substructure at the 2 TeV scale, and is competitive with the Tevatron despite 
the much lower c.m.s. energy 26 . Other candidates for new physics are possible 
new gauge bosons, e.g., a leptophobic Z'. Of course, high luminosity and preci¬ 
sion as well as a good knowledge of polarized and unpolarized parton densities 
and of the standard model ‘background’ are mandatory. For details, see 2,26 . 

5. Summary and Outlook 

With first data from RHIC hopefully starting to roll in soon, we can address 
many open, long-standing questions in spin physics like the longitudinally po¬ 
larized gluon density or transversity. With data from many different processes 
taken at high energies where perturbative QCD should be at work, a first 
global analysis of spin-dependent parton densities will be possible. At the end 
of RHIC we certainly have a much improved knowledge of the spin structure 
of the nucleon, and, perhaps, the next ‘spin surprise’ is just round the corner. 
Future projects like the EIC 27 , which is currently under scrutiny, would help 
to further deepen our understanding by probing aspects of spin physics not 
accessible in hadron-hadron collisions. The structure function g\ at small x or 
the spin content of circularly polarized photons are just two examples. 
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